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Neutron scattering from high-quality YBa 2 Cu 3 06.33 (YBC06.33) single crystals with a T c of 8.4 
K shows no evidence of a coexistence of superconductivity with long-range antiferromagnetic order 
at this very low, near-critical doping of p~0.055. However, we find short-range three dimensional 
spin correlations that develop at temperatures much higher than T c . Their intensity increases 
smoothly on cooling and shows no anomaly that might signify a Neel transition. The system 
remains subcritical with spins correlated over only one and a half unit cells normal to the planes. 

At low energies the short-range spin response is static on the microvolt scale. The excitations out 
of this ground state give rise to an overdamped spectrum with a relaxation rate of 3 meV. The 
transition to the superconducting state below T c has no effect on the spin correlations. The elastic 
interplanar spin response extends over a length that grows weakly but fails to diverge as doping is 
moved towards the superconducting critical point. Any antiferromagnetic critical point likely lies 
outside the superconducting dome. The observations suggest that conversion from Neel long-range 
order to a spin glass texture is a prerequisite to formation of paired superconducting charges. We 
show that while p c =0.052 is a critical doping for superconducting pairing, it is not for spin order. 


I. INTRODUCTION 

Despite numerous experimental and theoretical studies 
since the discover^ of high temperature superconductiv¬ 
ity (HTSC) in cuprates in 1986, there is still no consensus 
on the superconducting mechanism in these materials. 
With enormous advances in growth of high-quality sin¬ 
gle crystals, however, a qualitatively common phase dia¬ 
gram has emerged (Fig. [l]). The phase diagram exhibits 
different states of matter, many of which are revealed 
by neutron scattering 2 24 . Only when novel properties 
of HTSC cuprates are fully elucidated across the whole 
phase diagram will a full understanding of these mate¬ 
rials be gained. The search for new phases of matter is 
particularly important in the low doping region close to 
critical doping for superconductivity, p c , because it cre¬ 
ates the greatest challenge for theory. 

The undoped HTSC parent materials are believed to 
be Mott-insulators with static long-range three dimen¬ 
sional (3D) antiferromagnetic (AF) order below a Neel 
temperature T^r. The Cu 2+ spins are localized with 
strong AF exchange coupling within 2D Cu02 planes. 
The weak exchange coupling in the perpendicular direc¬ 
tion allows the 3D Neel order to occur at a finite Tjy. 
Initially, hole doping disrupts the network of correlated, 
localized Cu 2+ spins and suppresses antiferromagnetism. 
Further hole doping leads to superconductivity for dop¬ 
ing larger than the critical doping p c ^0.05 (underdoped 


region). The superconducting transition temperature T c 
increases with increasing doping and reaches a maximum 
at an optimal doping beyond which T c declines. Even 
though long-range AF order does not coexist with su¬ 
perconductivity in this region, strong spin fluctuations 
nonetheless survive well into the superconducting (SC) 
phase. Electronic properties more like conventional met¬ 
als are exhibited at very high doping (overdoped region) 
but also in the underdoped region by the observation of 
quantum oscillations and Fermi surface reconstruction 25 
in the range 0.08<p<0.15. Our results address a much 
lower doping region where superconductivity is weak and 
spin fluctuations are strong. 

The single-layer La 2 - cc Sr ;E Cu 04 (LSCO) cuprate fam¬ 
ily has been extensively studied by neutron scatter¬ 
ing 26 36 . The doped LSCO system intrinsically suffers 
from disorder very close to the Cu 02 planes in which 
both SC and AF occur. This is because doping involves 
substitution of La with Sr ions of different charge in lay¬ 
ers directly adjacent to the Cu02 planes. This chem¬ 
ical substitution inevitably introduces disorder into the 
Cu0 2 planes which in turn distorts the Cu02 planes. For 
the bilayer YBa2Cu306+cc (YBC06+x) system, where 
large and high-quality single crystals are available, dop¬ 
ing occurs by varying the oxygen content in the CuO 
chains which are located far from the Cu02 planes. The 
YBC06+X system thus provides an intrinsically less dis¬ 
ordered environment, where holes are introduced in the 
conducting plane by out-of-plane oxygen doping, that can 
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FIG. 1. Electronic phase diagram of YBC06+X is depicted as 
a function of hole doping. For the undoped parent material, 
there is an unambiguous transition to a long range 3D ordered 
AF phase at Tjy. For p > p c = 0.052 the system exhibits 
superconductivity below T c . How the long range AF order 
gives way to superconductivity by doping remains a central 
question in the physics of HTSC cuprate supercoductors. 


reveal how doping affects the planar magnetic and super¬ 
conducting properties. 

Muon spin rotation (fi SR) studies of underdoped 
YBC06+X have led to the claim^THU that a transition 
takes place to Neel order within the superconducting 
dome. As well, a transition to a frozen glass is claimed 
at an even lower temperature. One must emphasize 
that the fi SR technique is a local probe whereas neu¬ 
tron scattering allows one to determine the length and 
time scales of the magnetic structure as it evolves with 
temperature and doping. Since there is little consensus 
on the nature of the precursor phase from which the SC 
phase emerges, we have undertaken a systematic neu¬ 
tron scattering study of high-quality YBC06+X single 
crystals that lie close to the critical doping for super¬ 
conductivity (see Fig. [l]) . Our earlier experiments on 
YBC06.35 (T c =18 K, p=0.06) revealed that no long- 
range B rag g-ordered AF phase coexists with supercon¬ 
ductivity 22 24 [ The short-range static spin correlations 
observed for YBC06.35 leaves open the possibility that 
an AF quantum critical point or novel phase might occur 
at lower doping. We have therefore studied an even lower 
doped YBC06+X, namely YBC06.33 with T c of only 8.4 
K and doping p=0.055 near the critical p c =0.052. We 
will present clear evidence for the existence of spatially 
highly-correlated and slow spin fluctuations. Although 
the observations are qualitatively similar to higher doped 
YBC06.35 (T c =18 these are the first measure¬ 

ments on YBC06.33 with such a low T c , and we will 
show that even for such low doped yet superconducting 


materials, long-range antiferromagnetic order and super¬ 
conductivity do not coexist. 


II. EXPERIMENT 

The experiments were performed on large single crys¬ 
tals prepared at University of British Columbia (UBC) 
by a top-seeded melt growth technique 40 . The crystals 
each of volume ~1 cm 3 were well annealed and in the 
orthorhombic phase with room temperature lattice pa¬ 
rameters a=3.844 A, 6=3.870 A (for data analysis pur¬ 
poses we take a=b=(a + b )/2 =3.857 A), and c=11.791 
A. The superconducting transition temperature, T c , was 
determined to be 5 K after the crystal growth and ris¬ 
ing to 8.4 K after further annealing for several months 
at room temperature at Chalk River. The superconduct¬ 
ing transition temperature determined by AC Meissner 
susceptibility in a coil at 17 Hz remains very sharp for 
such a low doping (the transition width is only 2 K), see 
Fig.^a). T c is reduced to 2.5 K in a magnetic field of 7.5 
T || [110]. The sample hole doping p=0.055 was deter¬ 
mined from the observed T c in zero field and the empiri- 
cal formulcM 1-T c /T max =82.6[p-0.16] 2 with T max =94.3 
K which is also consistent with the c-axis lattice constant 
vs. doping relation observecP^ for YBC06+X. For the 
zero-field experiments, four crystals (4 cc in total with a 
mosaic spread of 2°) were co-aligned in the (HHL) plane. 
For the experiment in an applied magnetic field, one crys¬ 
tal (1 cc with a mosaic spread of 1°) was used with the 
field applied along [110]. 

Thermal neutron scattering experiments were made 
with the C5 spectrometer at NRU reactor at Chalk River 
Laboratories with pyrolytic graphite PG(002) crystals 
vertically focusing for the monochromator and flat for the 
analyzer with a fixed final energy of Ej=14.6 meV. Elas¬ 
tic and inelastic measurements of the sample aligned in 
the (HHL) plane were made in the temperature range 1.5 
K to 300 K, in magnetic field up to 7.5 T, up to 16 meV 
energy transfer in the (0.5 0.5 2) zone and up to about 40 
meV in the (0.5 0.5 5) zone. Fast neutrons were removed 
from the incident beam by a liquid nitrogen cooled sap¬ 
phire filter before the monochromator. Two PG filters 
with a total thickness of 10 cm in the scattered beam re¬ 
moved higher-order neutron wavelengths. The horizontal 
collimations were controlled with Soller slits as specified 
in the plots of data. A search revealed no oxygen su¬ 
perlattice peaks indicating that any oxygen ordering (i.e. 
Ortho-II and/or Ortho-Ill) in the CuO chains is short- 
ranged. Nonetheless a small a-b splitting was detected 
with x-rays, showing that the sample remains in the or¬ 
thorhombic phase. 

Cold neutron measurements were performed with the 
SPINS spectrometer at the NIST Center for Neutron Re¬ 
search (NCNR) and with PG(002) crystals vertically fo¬ 
cusing for the monochromator and flat for the analyzer 
with a fixed final energy of Ej=2.9 meV and a collima- 
tion of [guide 80' S 80' 80']. For energy transfers greater 
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FIG. 2. (a) The superconducting transition is observed at 

8.4 K at zero field and is reduced to 2.5 K at 7.5 T || [110]. 
The sample was located in a coil and the inductance of the 
coil, L, at 17 Hz was monitored as a function of tempera¬ 
ture using an inductance bridge. The inductance of the coil 
is related to the susceptibility of the sample, xs( T), through 
L = Lo(l + 47rxs(T)). When the sample becomes supercon¬ 
ducting, xs(T) becomes diamagnetic and hence the induc¬ 
tance of the coil will decrease, (b) Central mode intensity at 
E=0 at Qaf=(0.5 0.5 2) at zero field and at 7.5 T || [110] 
observed with thermal neutrons (resolution^ meV). Despite 
the change in T c , the central mode is not affected by the mag¬ 
netic field. Also shown is the intensity at zero field observed 
with cold neutrons (resolution^O.l meV) and with backscat- 
tering (resolution'll fieV). The intensity observed for each 
measurement is normalized to its maximum at low temper¬ 
atures. The temperature scale of the central mode intensity 
decreases when measured with better energy resolution. The 
solid line is a guide to the eye. 


than 2 meV a cooled Be filter was placed in the scattered 
beam. For energy transfers less than 2 meV a cooled Be 
filter was placed before the sample. Data was normalized 
from results at overlapping energies. 

To further investigate how the temperature depen¬ 
dence of magnetic correlations depends on the experi¬ 
mental energy resolution, we also carried out neutron 
scattering measurements with the NCNR’s high-flux 
backscattering spectrometer, HFBS, with an energy reso¬ 
lution of ^1 /ieV. Elastic scans were performed at a fixed 


final energy of Ej=2.08 meV using a Si(lll) monochro¬ 
mator mounted on a Doppler drive. Comparison of the 
temperature dependence of the observed magnetic elas¬ 
tic (central mode) data at (0.5 0.5 2), with thermal, cold 
and backscattered neutrons is shown in Fig. |2jb). 

The magnetic scattering function S(Q,lj), for momen¬ 
tum transfer Q, and energy transfer w, can be determined 
directly from the measured magnetic neutron scattering 
intensity by 

I(Q,v) oc f 2 (Q)B 2 (Q c )S(Q,u;) (1) 

where f(Q) is the anisotropic Cu 2+ form facto]®2 an d 
B(Q C ) = 2sin (^^z') is the independent bilayer struc¬ 
ture factor with Q c = 2ttL/c and z! m 1 — 2zc U 2 (see 
the Appendix). The model used for in the data 

analysis is discussed later in Section m The measured 
data is fitted to a function which is the convolution® of 
this model with the spectrometer resolution function to 
obtain the model parameters. All uncertainties and error 
bars are standard errors. 

The inelastic data are corrected for contamination of 
the incident beam monitor by higher wavelength harmon- 
icd®. A transverse acoustic phonon measured close to 
the (006) Bragg peak is used to put the observed in¬ 
tensities on an absolute scale, see Appendix. The total 
elastic magnetic moment is obtained by integrating over 
the observed elastic line in absolute units in energy and 
in momentum transfer along [HH0] and [00L] directions 
relative to the magnetic Brillouin zone centre. The re¬ 
sults are consistent with the magnetic moment derived 
from the analysis of several nuclear Bragg peaks. 

III. RESULTS 
A. Two energy scales 

The magnetic spectrum at the antiferromagnetic (AF) 
center Qaf=(0.5 0.5 2) at 1.5 K, measured with cold neu¬ 
trons with an energy resolution of 0.08 meV, is shown in 
Fig. [3] The average of the scattering observed at (0.3 0.3 
2) and (0.7 0.7 2) is used as background (BG 0 gp_p ea ^). 
This method of background subtraction is verified by 
constant-energy Q-scans at several energy transfers (see 
open circles Fig. [3] (b)). The spectrum exhibits two en¬ 
ergy scales: a very slow or static response characterized 
by a resolution limited elastic peak at zero energy (cen¬ 
tral mode) and a much weaker and broader spectral fea¬ 
ture. All spectra measured such as Fig.[3|b) indicate that 
there is no spin gap. At low energies the inelastic mag¬ 
netic scattering linearly decreases to zero with decreasing 
energy as may be expected from overdamped spin waves 
(paramagnons). The observed absence of dynamic spin 
spectral weight as energy goes to zero can only help for¬ 
mation of superconducting pairs for T<<T C . Any Neel 
order would instead have led to a growth in spin spectral 
weight as energy goes to zero, much as upon approach to 
a 3D Neel phase. 
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FIG. 3. Inelastic spectrum measured with cold neutrons at 
the AF center Qaf = (0.5 0.5 2) at 1.5 K. The average of the 
intensity at (0.3 0.3 2) and (0.7 0.7 2) is used as background 
(BG 0 ff_peak)- Scattering comprises an intense resolution lim¬ 
ited peak centered at zero energy transfer (a), and a much 
weaker broad component as seen in panel (b). The solid line 
is a fit to a resolution-limited Lorentzian at zero energy and to 
a broad relaxational Lorentzian with a ~3.4±0.3 meV relax¬ 
ation rate. The horizontal line in (b) is the energy resolution 
at zero energy transfer. The open circles in (b) are the peak 
height of the fitted constant energy H-scans along (0.5 0.5 
2), and they agree with two-point H=0.3 and 0.7 background 
subtraction. 


In the following we present the temperature and 
wavevector dependence of the central mode and the in¬ 
elastic feature. Our results demonstrates that the central 
mode represents only short-range magnetic correlations 
and that the dynamic overdamped fluctuations give rise 
to the broad inelastic feature. 


B. Quasielastic scattering: central mode 

Figs. [I] and [5| a) show typical elastic Q-scans through 
(0.5 0.5 2 ) along [HH0] and [ 00 L] directions observed at 
base temperature with a background subtracted at 80 K 
and 100 K respectively. The lack of temperature depen¬ 
dence above ^80 K justifies attribution of the growing 
scattering at low temperatures to the magnetic scatter¬ 
ing. This is confirmed by polarized neutron scattering®!. 


The scattering appears in form of a peak centered at 
the AF wavevector, Qaf =(0.5 0.5 L) with L=integer 
and broader than the resolution. To make a comparison 
of the observed peak widths with the resolution of the 
spectrometer, we directly measured the resolution by re¬ 
moving the PG filters and repeating the exact scans. The 
calculated resolution is in good agreement with the mea¬ 
surement (horizontal lines) as seen in Figs.[4jb) and (d). 
In addition, we made a comparison with the observed 
scattering from an undoped YBCO6.00 crystal with a 
Neel transition to a long range antiferromagnetic state 
at ~430 K. We find that the peak widths for the un¬ 
doped crystal are similar to the resolution and smaller 
than the ones in YBC06.33 along both directions. This 
comparison confirms in this low doped YBC06.33 su¬ 
perconductor the correlation lengths of the spins remain 
finite and long-ranged AF order is absent. The fact that 
the scattering is peaked at integer L-values, however, in¬ 
dicates the 3D spin correlations are parallel (ferromag¬ 
netic) between bilayers in adjacent cells along the c-axis. 
The spins in the two planes of a bilayer in one cell always 
remain antiparallel. 

The magnetic scattering function in Eq. |TJ 5 '(Q,cj), 
can be written as: 


S(Q,lo) = \F(Q)\ 2 g(cu) ( 2 ) 

where |F(Q )| 2 and g(uj) describe the momentum and 
spectral forms of the scattering, respectively. To extract 
the zero-frequency correlation lengths for both reciprocal 
space directions, we fitted the elastic data along [HH0] 
and [00L] to the convolution of the four-dimensional res¬ 
olution function 44 with the elastic magnetic scattering 
function of the form, 

5(Q,o;~0) = |F(Q)| 2 ^H (3) 


To avoid problems with the convolution, the delta func¬ 
tion was slightly broadened to 


5(iv) 


7 

uo 2 + q 2 


( 4 ) 


with 7 < A E res fixed where AE res is the energy resolu¬ 
tion, with 7 typically of order 1 fieV and 10 /ieV for cold 
and thermal neutron data sets, respectively. 

We first modeled the momentum dependent spectral 
form, |F(Q)| 2 , with a a set of independent Lorentzians, 


\m)\ = T 


Acm 


Aj 

l £2 


( 5 ) 


< h e ah 

where q ab = (H - 0.5, K - 0.5)(^), q C:i = (L — L t ){ 2 -f) 


with Li an integer ranging from -1 to 7, and ^ ab and £ c 
are correlation lengths in the ab-plane (assumed symmet¬ 
ric along the a and b in-plane directions) and along the 
c-axis, respectively. These Lorentzians describe the 3D 
correlations along L and their relative intensities, Ai, are 
allowed to vary independently. Because fitting of £ ab to 
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FIG. 4. Spin correlations from elastic scans along the [HHO] and [00L] through the AF center Qaf=(0.5 0.5 2). Data collected 
with cold neutrons (E/=2.9 meV) at 1.5 K are shown in panels (a) and (c) and with thermal neutrons (E/ = 14.6 meV) at 3 K 
in panels (b) and (d). The horizontal lines are the calculated resolution Bragg widths (FWHM). To make a direct comparison 
with the resolution of the spectrometer at this momentum transfer, we removed the PG filters and repeated the scans (red 
hexagons). In addition, we measured the magnetic scattering from an AF Bragg ordered single crystal of YBCO6.00 with a 
transition to a long-range AF order (green stars). The magnetic scattering for YBC06.33 is wider than both the measured 
resolution of the spectrometer and the magnetic Bragg peak observed for YBCO6.00. Gaussian FWHM for each scan is given 
in brackets in the panels b and d. This indicates that only short-ranged magnetic correlations exist in YBC06.33. The fits are 
obtained by convolution of the 4D resolution with an elastic Lorentzian function with isotropic q-width in the a*-b* plane and 
a larger width along c*. 


an H-scan involves resolution convolution over a known 
L width (f c ), and vice versa, we iterated to ensure that 
a consistent set of and £ c were fitted to the two data 
scans (Figs. [4] and §a)). We extract the intrinsic low 
temperature elastic spin correlation lengths at 1.5 K to 
be £ a 6=110±20 A in the ab-plane and £ c =17±l A for 
all peaks along the c-direction. Both cold and thermal 
data gave the same values within the error bars after 
the resolution effects are included. Directly from the ob¬ 
served widths of the data along [HHO] and along [00L] 
directions a correlation length of about 40-50 A is ob¬ 
tained for in-plane correlations (l/£ a & = v/2AH(27r/a)) 
and about 11-17 A for correlation length along the c- 
direction (l/£ c = AL(27r/c)), where AH and AL are the 
observed Lorentzian HWHM in r.l.u. along H and L, 
respectively. 

We may also extract the c-axis correlation length, £ c , 
by assuming that the bilayer amplitude decays exponen¬ 


tially along the c-axis as exp(— r z /£ c ). The Fourier trans¬ 
form of this correlation is exactly giveiP^ by 

iFfO'll 2 = sinhg~ 1 c/2) cosh(C~ 1 c/2) 

1 + llbtlh sinhg^c/ 2 ) 2 + sin(Q c c/2) 2 

For large £ c , Eq. [6] gives a set of Lorentzians along L 
centered at each integral Li but with a single £ c (of 
the form given by Eq. |5j. This coupled bilayer for¬ 
mula gives the intensity of the peaks along L through 
the bilayer structure factor B 2 (Q C ) and eq. [l] The £ c 
=16 zb 1 A extracted from a fit of the coupled bilayer 
model agrees with the value obtained using the indepen¬ 
dent Lorentzians of Eq. (5). The coupled bilayer model 
used in Ref. 22} is the square of Eq. [6] and hence the c- 
axis correlation lengths reported there correspond to a 
Lorentzian-squared model. 
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FIG. 5. (a) Approach to 3D via out-of-plane elastic spin cor¬ 
relations along [00L] through [0.5 0.5 L] at 3 K where data at 
100 K is used as background. Fits are shown of resolution- 
convoluted Lorentzians centered at integer L (solid red line) 
and to the coupled bilayer model (dashed green line). For the 
convolution with resolution the in-plane correlation lengths 
are taken from elastic scans along [H H 0] through (0.5 0.5 
2), similar to the data shown in Fig. [4] (b) The L-scan at 
a finite E=4.14 meV energy transfer shows that spin exci¬ 
tations mainly originate from independent bilayers (dashed 
green line). The data is corrected for background taken as 
the average of the observed intensity at (0.3 0.3 L) and (0.7 
0.7 L), BG 0 fpp ea k-. (c) A simulation of scattering along the In¬ 
direction for uncoupled bilayers (dashed green line) and ferro- 
magnetically coupled bilayers (solid red line). Horizontal lines 
at L=2 in (a) and (b) give the calculated resolution widths 
in L (FWHM) at E=0 and 4.14 meV. In (a) and (c) the in¬ 
cipient 3D pattern is larger than in (b) because the static 
susceptibility is much larger as the energy approaches zero. 


Fig. [^a) shows that the coupled bilayer model de¬ 
scribes the data almost as well as the independent- 
Lorentzian model but underestimates the intensity ratio 


of the peak at L=1 to that at L=2. Since the inten¬ 
sity recovered somewhat when cold neutrons were used, 
we attribute the difference to unknown resolution effects. 
The fit is equally good to a Lorentzian squared with a cor¬ 
relation length along the c-axis of ^7.8±0.5 A, shorter 
as expected. 

In Fig. [6] we show the properties of the central mode in 
the L-E and L-H planes. From Gaussian fits we plot in 
panels (b) and (d) the peak height, and the energy and 
momentum widths. As L^2, the central mode amplitude 
grows while the relaxation rate slows to the resolution 
limit of 0.08 meV. Right at the integer-L momentum, the 
spins are frozen on the neutron timescale. The smooth 
behaviour in peak height, energy and momentum width 
as L tends to the 3D point L=2 shows that there is no 
Bragg feature hidden within the central mode correla¬ 
tions. Thus the spin correlations are short range and no 
Neel order is present. The near-critical scattering does 
not follow a Lorentzian form, for the width in H would 
then continue to grow as L departed from L=2. Instead 
the H width becomes constant far from L=2. This be¬ 
haviour shows that the quasielastic scattering far from 
the L-integral peaks adopts a rod-like form typical of 2D 
systems with a constant width in H. 

Figs.[7][a,b,c) show the temperature dependence of the 
quasielastic scattering near (0.5 0.5 2) by scans in energy, 
H and L. The fit parameters including central mode am¬ 
plitude, the quasielastic energy width FWHM obtained 
from a gaussian fit, and the inverse correlation lengths 
obtained from the convolution of resolution with Eq. [5] 
are shown in Figs. [7|[d,e,f). We find that the quasielastic 
amplitude grows monotonically on cooling before satu¬ 
rating below about 15 K. There is no anomaly that might 
signal a transition to long-range AF order. In addition, 
the transition to the SC state at T c =8.4 K does not affect 
the temperature dependence of the central mode indicat¬ 
ing that superconductivity and AF spin coupling develop 
independently. 

The inverse correlation lengths along both directions, 
£ a 5 and f c , remain finite at all temperatures. The ^ 
gradually decrease on cooling while £ c exhibits a more 
pronounced decrease below ^30 K. The spins are at¬ 
tempting to find a state similar to the AF state below 
30 K but are prevented by hole doping. Both inverse 
correlation lengths may hint at the presence of a mini¬ 
mum around ~15 K. As well the central mode dynamics 
slows on cooling until reaching the resolution limit. A 
similar narrowing of the dynamic width was observed 22 
in YBC06.35 (T c =18 K) but the increase of the spin 
correlation lengths below 15 K on cooling in YBC06.33 
(T c =8.4 K) is a new effect (Figs. [7je,f)). The growth 
of the central mode amplitude cannot be accounted for 
by the decrease in its extent in momentum-energy space, 
£ab Cr 1 AE 1 , all of which decrease from 80 K to 15 K. 

A comparison of the temperature dependence of the 
central mode intensity measured with different resolu¬ 
tions shown in Fig. [2jb) indicates that the temperature 
scale of the central mode depends on the energy resolu- 
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FIG. 6. (a) Quasielastic energy scans of the central mode as a function of L close to AF center Qaf=(0.5 0.5 2) at 1.5 K. (b) 
Central mode peak height and energy FWHM obtained from scans in (a) by fitting to a Gaussian function, (c) Elastic H-scans 
of the central mode as a function of L at 3 K. (d) Central mode peak height and FWHM in H obtained from (HH2) scans in 
(c). In (b) and (d) dashed lines are the calculated resolution and the solid lines are guides to the eye. 


tion at which the measurement is performed. When the 
central mode is probed with high-resolution backscatter- 
ing with an energy resolution of only ^1 /ieV, the scat¬ 
tering appears on cooling only below ^40 K compared 
to ^80 K for thermal neutrons with broader energy res¬ 
olution of ^1 meV. Hence all thermal and cold neutron 
measurements are extremely resolution limited, e.g., res¬ 
olution of 0.08 meV is 80 times wider than the maximum 
intrinsic width of 1 /ieV of the central mode. The depen¬ 
dence of the central mode intensity vs. T on resolution is 
similar to the behaviour observed in spin glasses at low 
temperature^^!!. 


C. Inelastic spectra 

Typical inelastic spectra measured by const ant-Q en¬ 
ergy scans at (0.5 0.5 2) are shown in Fig. J8^a) at dif¬ 


ferent temperatures. The spectra are corrected for the 
wavelength feedthrough in the monitor 45 . The average of 
the scattering at (0.3 0.3 2) and (0.7 0.7 2), BG 0 gp_p ea ^, 
is used for background subtraction. This method of 
background subtraction is confirmed by constant-energy 
(HH2) scans through the AF center Qaf=(0.5 0.5 2) 
resulting in the same peak intensity (open circles in 
Fig. |8Fa)). We have used the phonon calibration method 
(see Appendix) to put the observed scattering on an ab¬ 
solute scale (right-hand axis in Fig. |8ja)). 

The background corrected data is fitted to the convo¬ 
lution of the 4D resolution function with the model of 
Eq.| where g{uj) is given by 


9 M 


7 Acm_ 

w 2 + q 2 


+ [1 + n{uj)\ 


ujTA 


ex 


w 2 + r 2 


( 7 ) 


where [1 + n(w)]=l/[l — exp(—is the Bose 
population factor. The first term describes an energy 
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FIG. 7. The temperature dependence of the intensity around Qaf =(0.5 0.5 2) in (a) energy scans around zero energy 
transfer, (b) in elastic (HHO) scans and (c) in elastic (00L) scans through (0.5 0.5 2). The data is corrected for background by 
subtracting the data at 80 K in (a) and (b) and at 100 K in (c). The dotted horizontal line in (d) is the calculated resolution 
width (FWHM). No signature of a Neel transition down to 1.5 K is found in the temperature dependence of the central mode 
amplitude. The central peak is consistent with being resolution limited in energy at low temperatures. As seen in panels (e) 
and (f), the spatial correlation lengths in the ab-plane and along the c-direction lengthen gradually on cooling but remain finite 
at the lowest temperatures. 


resolution-limited central peak forced to be elastic by 
setting 7 A E res . The second term represents spin 
diffusion described by a modified Lorentzian defined by 
its relaxation rate, T. The concurrent Lorentzian mo¬ 
mentum dependence (Eq. [5| displays a memory loss over 


a time ( 1 /T) and over long but finite corresponding cor¬ 
relation lengths in the ab-plane and very short distances 
along the c direction. The solid lines in Fig. | 8 ja) show 
that the model can fit the data reasonably up to 10-15 
meV. We find that at higher energy transfers the data 
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FIG. 8. (a) Temperature dependence of inelastic spectra observed at the AF center Qaf= (0.5 0.5 2). At each temperature the 

averaged scattering measured at (0.3 0.3 2) and (0.7 0.7 2), BG 0 fj_p ea k, is used for the background subtraction (solid circles). 
This method of background subtraction is confirmed by const ant-energy (HH0) scans through the AF center Qaf = (0.5 0.5 
2) resulting in the same intensity (open circles). The solid lines are fits to a broad damped response with a relaxation rate, 
r, given by Eq. [7] Note that for these fits, the data within the elastic resolution region is removed, (b) The temperature 
dependence of the model parameters Xo=A ex and F. A Curie-Weiss high temperature trend for F is shown by a dashed line, 
and that for the static susceptibility by a solid line. We find that a quadratic temperature dependence describes the relaxation 
rate, T(T)=a + bT 2 , over the entire temperature range (the dash-dotted line). 


falls less rapidly than the model. Extra scattering at en¬ 
ergies larger than 15 meV could form a precursor to the 
almost constant local susceptibility expected for 2D spin 
waves at larger energies. At larger energies, E^20 meV, 
phonon contributions are known 18 ! to occur. Note that 
the width in q is constant for E<10 meV (Fig. [Tljb)) so 
that the spectrum at Qaf=(0.5 0.5 2) is proportional to 
the local susceptibility. 

The temperature dependence of the model parameters 
are shown in Fig. [8^b). The amplitude of the modi¬ 
fied Lorentzian Xo=A ex is proportional to the suscep¬ 
tibility of the spin fluctuations. This dynamic suscep¬ 
tibility gradually increases on cooling until it saturates 
and turns over at low temperatures. The remainder 
of the spin spectrum lies unresolved within the central 
mode discussed above. The data, however, shows that 
the dynamic AF correlations are still present within the 
bilayer at temperatures much higher than the temper¬ 
ature scale for appearance of the static central mode. 
The solid line in Fig. [8jb) is the fit to a Curie-Weiss 
behaviour, %o = C/(T + &c), with Oc ^32±8 K. As 
seen the data deviates from this behaviour below ^30 K. 
The phenomenological Curie-Weiss temperature is con¬ 


sidered the temperature scale for the appearance of short- 
ranged AF correlations at high temperatures in the para¬ 
magnetic state of geometrically frustrated magnets 51 and 
some heavy fermiond^. 

Fig. gb) shows the relaxation rate decreases al¬ 
most linearly on cooling down to about ^80-100 K 
(dashed line) and eventually saturates to a constant 
value below ^30 K. A linear temperature dependence 
for the relaxation rate is common in geometrically frus¬ 
trated system^ 51 ! and it is also observed for lightly 
doped La 2 Cuo. 94 Lio.O 6 and Lai. 95 Bao.o 5 Cu 04 54 . In 
YBC06.33 we find instead that a quadratic tempera¬ 
ture dependence describes the relaxation rate, T(T)=a 
+ bT 2 , over the entire temperature range as shown with 
the dash-dotted line in Fig. [8|b). 

Our analysis indicates that the inelastic spectrum can¬ 
not be described by a linear spin wave spectrum convo¬ 
luted with the 4D resolution. In particular the inelas¬ 
tic peak cannot arise from the known vertical resolution. 
As shown later in Fig. [Tl] and discussed in relation to 
Fig. [l3j the finite extent in momentum of the spin cor¬ 
relations prevents divergence of the low energy suscepti¬ 
bility. Right at the staggered wavevector this results in a 
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FIG. 9. A comparison at (0.5 0.5 2) of the intensity of the cen¬ 
tral mode with that of the low energy fluctuations at E=0.5 
meV as a function of temperature. A spectral transfer from 
dynamic fluctuations to elastic correlations is observed as the 
system is cooled down to its ground state. Solid lines are 
guides to the eye. The characteristic temperature of ~30 K 
for the depopulation of the fluctuations corresponds to the 
relaxation rate of meV and not to the 0.5 meV fluctua¬ 
tion energy. The loss of fluctuation strength appears in the 
growth of the ground state represented by the central mode. 



FIG. 10. Temperature dependence of the integrated intensity 
of elastic and inelastic scattering. As the inelastic scattering 
decreases on cooling, the elastic component increases. Solid 
and dashed lines are guides to the eye. 


spin suppression leading to an overdamped peak centered 
on ^3 meV. 

A direct comparison of the temperature dependence 
of the intensity of low energy fluctuations with that of 
quasielastic scattering is shown in Fig. [9] where the in¬ 
tensity of the central mode and the excitations at 0.5 


meV are plotted vs. temperature. Both data sets are 
collected with cold neutrons and are corrected for back¬ 
ground. The central mode peak height grows on cooling 
while the low-energy fluctuations diminish. This shows 
that the excitations are indeed thermal fluctuations of 
the quasistatic ground state at low temperatures. This 
correlated behaviour also suggests that they arise from a 
single spin phase linked together by the hole doping, yet 
independent of the charge pairing. 

Fig. [lO] shows the integrated intensities of the inelas¬ 
tic and elastic scattering from thermal neutron data as 
a function of temperature. The inelastic integration is 
obtained from the constant-Q energy scans by numerical 
calculation of the area underneath the observed scatter¬ 
ing from 1 meV up to 16 meV. This is justified since 
the Q-scans at different energies in this low energy range 
exhibit similar widths (see Fig. 11). The elastic integra¬ 
tion is obtained from elastic Q-scans in both H- and In¬ 
directions simply by multiplying the observed amplitude 
by the observed widths in both in-plane directions [HH0] 
and [-HH0], and by the width in the [00L] direction. The 
inelastic scattering decreases on cooling mirroring the in¬ 
crease in the elastic scattering. In this figure we have 
also shown the sum of the integrated intensities of elastic 
and inelastic components. The sum is temperature inde¬ 
pendent below ^100 K as expected from conservation of 
the total moment. The decrease at higher temperatures 
arises because the scattering extends to much higher en¬ 
ergies than our upper limit of integration of 16 meV. 

We have also performed constant-energy Q-scans along 
[HH0] in the energy range from zero to ^40 meV and 
along [00L] from zero to 4 meV at different temperatures 
as shown in Figs. [TTl a), [I2| and [5|b). The half-width- 
at-half-max (HWHM) in A^ 1 , obtained from fitting the 
measured Q-scans along [HH0] to the convolution of res¬ 
olution with the Lorentzian form of Eq. [5j is shown in 
Fig. [lTJb) as a function of energy. For any nonzero en¬ 
ergy where the central mode does not contribute, the 
fluctuations are three times as broad in H as the static 
response. Their widths remain constant at ^0.03 A -1 up 
to 15 meV above which they increase and tend to follow 
the spin-wave velocity measured 43 for the AF long-range 
ordered insulator YBC06.15, co=600 meVA. 

Constant energy Q-scans along [0.5 0.5 L] shown in 
Figs. §b) and [12] indicate that the integer-L centered 
coupled bilayer behaviour observed at zero energy is 
greatly reduced as energy increases. By 2 meV the dy¬ 
namic correlations have lost the three-dimensional char¬ 
acter leaving only the pattern for independent bilayers 
centered on L=1.7. 

In order to investigate whether the relaxation rate rises 
with Q-Qaf, we measured inelastic scattering for nearby 
wavevector (0.52 0.52 2). Again we took the averaged in¬ 
tensity at (0.3 0.3 2 or 5.1) and (0.7 0.7 2 or 5.1) as 
background, BG Q gp_p ea ^. Fig. 13 shows that the weaker 
off-peak fluctuations display the same relaxation rate in¬ 
dependent of in-plane momentum H. Thus there is no 
dispersion associated with the low energy fluctuations in 
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FIG. 11. (a) The H-dependence of the scattering at different 
energy transfers around (0.5 0.5 2) at 1.5 K measured with 
thermal neutrons, (b) The resolution corrected HWHM in 
H of the inelastic scattering is constant and independent of 
temperature and energy from 1 to 15 meV and three times 
the elastic width (0.01 A -1 ). The slope of the inverse spatial 
correlation length at high energies is similar to the spin-wave 
velocit}^ (broken line) of insulating YBC06.15. 


YBC06.33. However the peak intensity at H=0.52 has 
fallen by only 1.3, much less that than the factor ^6 ex¬ 
pected for a Lorentzian dependence in q compared with 
the chosen offset of dH=0.02 r.I.u equivalent to dq =0.045 
A -1 where dqx £ a 5 = 4.5. The fact that the dynamic 
susceptibility holds up as q is displaced from Qaf shows 
the tendency, even at this very low doping, for the spin 
strength to be slightly weighted towards incommensurate 
wavevectors as also suggested by Fig. [Tljb). Moreover, 
for all energy scans up to 40 meV there is no sign of 
any other spectral feature that can be associated with 
the so-called resonance observed in samples with higher 
doping. 



FIG. 12. The L-dependence of the scattering at 1.5 K at 
different energy transfers measured with cold neutrons. The 
coupled bilayer modulation declines with energy and the L- 
dependence of the scattering becomes close to the indepen¬ 
dent bilayer form for E>2 meV. 


D. Magnetic field dependence 

Since one could argue that the transition to the super¬ 
conducting state prevents further growth of the central 
mode at low temperatures (i.e. prevents a transition to a 
long-range AF state), we suppressed the superconductiv¬ 
ity by applying a magnetic field of 7.5 T parallel to the 
[110] direction. T c is reduced to 2.5 K from 8.4 K at 7.5 T 
as shown in Fig. [2^a) . Detailed measurements of the elas¬ 
tic scattering along [HH0] and [00L] directions through 
(0.5 0.5 2) at 1.5 K are shown in Fig.[l4|and the temper¬ 
ature dependence of the central mode intensity is shown 
in Fig. [2](b) . Our standard method of high-temperature 
background subtraction was used. These measurements 
demonstrate that the central mode (intensity, q-profile 
and peak widths) and its temperature dependence are un¬ 
affected by the application of 7.5 T magnetic field, even 
though the system remains in its normal phase down to 
2.5 K. It appears that the static spins in the central mode 
evolve separately from the charges that pair below 8.4 K. 


In Fig. 15 [a-e) we show that the spin fluctuations both 
in q and in E are also unchanged by the applied field of 
7.5 T. These observations and the absence of an anomaly 
in both elastic and inelastic scattering at T c confirm 
that superconducting properties play an unobservable in¬ 
fluence on the very strong spin glassy correlations and 
their excitations. We do not probe the region near 0.8 
meV, where a field enhancement at 6 T was seen 5 ^ in 
YBC06.35. This region is close to the strong elastic line 
and is hard to probe with the lower-resolution thermal 
neutrons. There is possibly very tentative evidence of a 
field enhancement near 1 meV (Fig. [l5fc). 

The lack of field dependence in both elastic and in¬ 
elastic magnetic scattering in YBC06.33 contrasts with 
the LSCO systenP^^l In LSCO, for underdoped sam- 
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FIG. 13. Inelastic spectrum measured at (a) AF center 
Qaf=(0.5 0.5 2) and (b) off the AF peak position (0.52 0.52 
2) by twice the inverse correlation length (see Fig. [7|. In 
both panels high energy data is collected at (0.5 0.5 5.1) and 
corrected for the form factor. Again the average intensity ob¬ 
served at (0.3 0.3 2 or 5.1) and (0.7 0.7 2 or 5.1) is used as 
background, BG 0 fj_p ea k. Modified Lorentzian with the same 
T can describe both sets of data. Only the intensity is reduced 
at (0.52 0.52 2) compared with (0.5 0.5 2). 


pies with a Sr concentration of x<0.1 5, th e incommensu¬ 
rate magnetic scattering is enhanced 56 57 by application 
of a magnetic field. In the low doped superconducting 
La 2 - cc Ba cc Cu 04 with x=0.095 it was founcP^that the in¬ 
commensurate magnetic scattering is unaffected by mag¬ 
netic fields up to 7 T applied along the c-axis. The lack 
of a field effect in this system compared to the presence of 
strong field dependence in LSCO was used to suggest— 
that the field effect is not a universal property of cuprate 
superconductors. 

Experiments on more highly-doped YBC06.45 with 
T c =35 K have revealed 1 ^ a strong enhancement of static 
incommensurate magnetic order at low temperatures by 
a magnetic field of ^15 T applied mainly perpendicular 
to the Cu0 2 planes. It was also found that the field re¬ 
duces the amplitude of the inelastic response around ^3 
meV compensating for the spectral weight accumulated 
in the elastic peak. A similar field-induced suppression of 
the intensity is well known 65 for the resonance mode in 




FIG. 14. The elastic scattering does not change in an applied 
field of 7.5 T parallel to [110] as seen in scans (a) along [HH0] 
and (b) along [00L] through (0.5 0.5 2). The horizontal lines 
are the calculated resolution Bragg widths (FWHM). 


more highly doped YBC06.6. Based on these results 
it has been suggested 15 ! that the field enhanced mag¬ 
netic superstructure is expected to drive a reconstruc¬ 
tion of the Fermi surface which ultimately could explain 
the unusual Fermi surface topology revealed by recent 
quantum-oscillation experiments^. Recent experiments 
on lower doped YBC06.35 have instead shown 55; a mag¬ 
netic field induced enhancement of the low en ergy spin 
fluctuations. This enhancement was suggested^ 5 to be 
the result of free spins located close to the hole-rich re¬ 
gions. Clearly further experiments are required in order 
to determine whether there is a common magnetic field 
dependence for the static and dynamic scattering as a 
function of doping in YBCO and other cuprates, but the 
initial results suggest otherwise. 


IV. DISCUSSION AND CONCLUSIONS 

We find that in YBC06.33 only short-range AF corre¬ 
lations are present. Despite the near-critical doping, the 
system lies inside the superconducting phase but outside 
any antiferromagnetic ordered phase. The spins fluctu¬ 
ate on two energy scales: one a damped spin response 
with a ^3 meV relaxation rate and the other a resolution 
limited and intense peak at zero energy transfer (central 
mode). Even though we observe no transition to a long- 
range AF state, we demonstrate that the spins are highly 
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FIG. 15. (a)-(d) Constant-energy Q-scans at 1.5 K in an applied field of 7.5 T parallel to [110] direction compared to zero-field, 
(e) The inelastic spectrum at 1.5 K measured at (0.5 0.5 2) in an applied field of 7.5 T field parallel to [110] compared to zero 
field. Background again was taken as the average intensity at (0.3 0.3 2) and (0.7 0.7 2), BG 0 g-_p ea k-. No held dependence is 
observed for inelastic scattering similar to elastic scattering shown in Fig. [14] 


correlated and exhibit an imprint of three-dimensional 
ordering with finite correlation lengths. 

The central mode develops gradually on cooling but 
does not diverge and the system remains subcritical. We 
find no evidence for a phase transition such as a Neel 
anomaly (Figs. [2jb) and [7]). The intensity of the cen¬ 
tral mode saturates on cooling below ^10-15 K, with a 
concomitant minimum in the central mode peak widths. 
The presence of zero energy spin correlations over a fi¬ 
nite wavevector fraction of the Brillouin zone is enough 
to prevent a transition to 3D long-range order. 

The central mode is centered on commensurate AF po¬ 
sitions (0.5 0.5 L) with integral L indicating the develop¬ 
ment of weak 3D spin correlations along the c direction. 
The AF elastic correlations extend over only ^30 cells 
in the ab-plane and ^1 cell along the c-axis. The short- 
range nature of the correlations can be understood by 
the AF frustration due to the doped holes. The doped 
holes produce extensive regions which break up the AF 
coupling and likely create ferromagnetic correlations be¬ 
tween spins connected by oxygen neighbours in the Cu 02 
planes where holes reside. This frustration leads to a 
spin-glass state. Such frustration becomes more signifi¬ 


cant at low temperatures when the ferromagnetic bonds 
become frozen and thereby exert a nonzero average field 
on the regions of AF correlated spins. The correlations 
are broadened on warming and the integral-in-L pattern 
becomes centered on the maxima of the bilayer structure 
factor, so indicating that 3D coupling between bilayers 
is vanishing. The avoidance of a quantum critical tran¬ 
sition to 3D antiferromagnetism (divergence of the cen¬ 
tral mode intensity) is possibly related to breaking of the 
weakest coupling in the system, the antiparallel orienta¬ 
tion of the spin of the upper member of the bilayer in 
one cell with the lower member of the bilayer in the next 
cell up along c. This is also evident from the shorter 
correlation length along c than the ab-plane. 

There have been several suggestions for the local tex¬ 
ture that gives rise to the central mode. Haas et al. sug¬ 
gested^ skyrmions but we find the skirt around the base 
of the central peak is hard to differentiate from our excel¬ 
lent Lorentzian fits. Some have suggested^ 72 spin spi¬ 
rals. Aharony et al. suggested 73 spin canting in the plane 
as frustration that destroys the Neel order at such a low 
doping. It is known that a doped hole, forms a Zhang- 
Rice singled that neutralizes a Cu spin. Since the doped 
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TABLE I. Correlation lengths vs. doping for YBC06.33 
compared to YBC06.35. Data for YBC06.35 is taken from 
RefP. 


Material 

doping 

Tc (K) 

Cab (A) 

£c(A) 

YBC06.33 

0.055 

8.4 

110±20 A 

17±1 A 

YBC06.35 

0.06 

18 

42±5 A 

8±2 A 


hole extends over four or more oxygen sites, it is far more 
effective in destroying the long-range AF, reducing Tjy 
and magnetic moment, than simply a magnetic-site di¬ 
lution. To date the local texture around a hole has not 
been established. 

The in-plane correlation length of about 100 A is much 
larger than the hole spacing expected for a 2D system 
with doping p=0.055, l=a/y/p= 16 A. The correlation 
seems to lie closer with p=0.055 to a critical behaviour 
l— a/yjp — p c =7 0 A. This is not unexpected as the spin 
system is becoming more three-dimensional as p c is ap¬ 
proached. The trend with doping of the correlation range 
is shown in Fig.[l6j The inverse correlation lengths along 
both directions track T c without any sign of impending 
divergence at p c . It is remarkable that reducing dop¬ 
ing by only about 10% more than doubles the correla¬ 
tion lengths (see Table. |T| indicating high sensitivity of 
the static magnetic correlations to doping in this region. 
The closely coupled in-plane and out-of plane spin evolu¬ 
tion we find suggests a highly organized structure rather 
than doping disorder. We also note that the develop¬ 
ment in q, cj, and T of spin fluctuations is reminiscent 
of approach to a 3D phase transition and unlike what 
may occur in a system with disorder or random impuri¬ 
ties. From a linear extrapolation of the l/£ c observed for 
YBC06.35 and YBC06.33 (Fig. [l6|, we suggest a true 
3D long-range AF order to take place outside of SC dome 
at a doping lower than p c =0.052. This suggestion is lim¬ 
ited by the experimental accuracy, and is only possible 
because of the superb control of doping in the growth 
procedure used by the UBC group. Moreover there is 
little sign of a dynamic slowing of the spins as doping is 
reduced from YBC06.35 to YBC06.33. Evidently p c is 
a critical doping for the SC pairing, but not for the AF 
spin order. 

The quasielastic moment, being the integral in q of the 
central mode at base temperature, declines rapidly with 
doping. Our results shown in Fig.[l7]lie on a single trend 
line with doping that is the same for the YBCO and 
LSCCP systems. The value of the observed moment 
is derived by calibrating magnetic intensities from the 
observed intensities for several nuclear Bragg peaks 23 75 . 
We have verified that integration of the elastic line when 
put on absolute scale by phonon calibration also results 
in a similar value to within error bars (0.08T0.02 /ib 
compared to 0.11T0.02 ps from Bragg peaks). Their 
average is plotted in Fig. [17| 

Using high resolution probes (cold neutrons and the 
backscattering method) we have shown that while the 
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FIG. 16. The trend with doping of the correlation lengths in 
the ab-plane and along the c-axis are shown (right-hand y- 
axis). Doping dependence of superconducting transition tem¬ 
perature is also shown (left-hand y-axis). Data for YBC06.35 
and YBC06.5 is taken from Ref. 1231 As seen the correlation 
lengths become smaller as doping increases into the supercon¬ 
ducting dome. For a small transition temperature of only 8.4 
K in YBC06.33, the correlation lengths remain finite. 
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FIG. 17. The quasielastic moment, being the integral in q 
of the central mode, declines rapidly with doping. Data for 
LSCO system is taken from Ref. Data for YBC06+X with 
p > 0.06 is taken from Ref. M 


central mode is never resolution limited in q (i.e. is not 
Bragg-like), it becomes resolution limited (less than ~1 
/ieV) in energy domain at low temperatures. Thus the 
observed magnetic scattering appears static on a time 
scale of order of ^10 -9 s. The central mode, being es¬ 
sentially elastic yet extending over a finite momentum 
range, breaks the relation of energy to momentum via 
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velocity used to model highly-doped metallic cuprate!®. 
The short-range correlations within and between planes 
reveal a ground state of frozen, sub-critical, 3D-enhanced 
spin correlations. 

We find the that the temperature scale for the growth 
of the correlations depends on the energy resolution of 
the probe (see Fig. |2jb)). This behaviour observed for 
more conventional spin-glass systems'^® an d a i so seen 
in the underdoped spin-glass phase of LSCO 32 ^ is taken 
as an indication that the observed scattering is not truly 
elastic. This is because a probe with a tighter energy res¬ 
olution is more sensitive to slower fluctuations. It cannot 
detect the scattering at higher temperatures since the 
fluctuations are then faster than the time scale set by 
the energy resolution. The quasielastic characteristic of 
the scattering therefore can explain why the temperature 
scale of the central mode intensity depends on the energy 
window of the measurement. In YBC06.33 magnetic 
correlations fluctuate with multi-relaxational rates and 
avoid a transition to a static order. From this behaviour 
and the fact that the correlation lengths remain finite, 
we deduce that the magnetic ground state of YBC06.33 
is spin-glass like. The spin-glass behaviour coexists with 
superconductivity and is not affected by the transition 
to the SC paired state b elow T c nor by magnetic fields 
in contrast to the LSCO®® and LBCCP 7 ^ systems. It 
appears that ordering of spins and SC paired charges 
proceed separately. 

We do not observe a difference in the central mode in¬ 
tensity between field cooled and zero-field-cooled proto¬ 
cols as found in the spin-glass PbFe 1 / 2 Nb 1 / 2 03 by Chillal 
et af 78 . This may be because of the much larger exchange 
in YBC06.33 so that field effects are reduced below the 
resolution limit. Nonetheless our saturation of the peak 
intensity (below ~ 10-15 K) does mirror the spin-glass 
behaviour observed by Chillal et al. Hence we view the 
spin-glass freezing temperature of YBC06.33 as equal 
to ~10-15 K since below this temperature, the intensity 
and the peak widths saturate to a low temperature value. 
The spin-glass freezing here is different from that in ge¬ 
ometrically frustrated spin glasses but the concept can 
be applied to low doped YBC06.33 since the spins are 
frozen in time and have no long-ranged order. 

The absence of a transition to a long-range AF state 
differs from /i SR result^® 39 in underdoped YBC06+X, 
Ca-doped YBCO and LSCO. The phase diagram con¬ 
structed from such /i SR measurements includes two mag¬ 
netic transitions within the superconducting dome (a 
Neel transition at T^ and a transition to a frozen spin 
glass at a lower temperature), where an ordered AF state 
is suggested to coexist with SC below Even though 
the sensitivity of fi SR to slow spin fluctuations allows one 
to determine the onset temperature of magnetic correla¬ 
tions more easily and accurately, it only provides infor¬ 
mation on magnetic correlations integrated over q-space 
and therefore other complementary techniques such as 
neutron scattering are required to identify their spatial 
extent. We not only observe a gradual growth of the AF 


signal on cooling but also find that the correlation lengths 
in all directions remain finite to the lowest temperatures 
measured as is evident from a comparison of the observed 
magnetic peak widths with the experimental resolution 
and from a sample with long-range order (Fig. [4|. We 
note that the gradual onset temperature for appearance 
on the microvolt scale of the quasielastic central peak 
we observed with backscattering is comparable with T n 
determined from /iSR measurements where the onset is 
sharp in T. 

We have determined that the magnetic correlations 
occur in form of broad but commensurate peaks. The 
LSCO family formd^MHS] a spin-glass phase with short- 
range but incommensurate AF correlations coexisting 
with the SC phase for doping levels close to p c . The in¬ 
commensurability in the spin-glass phase scales linearly 
with doping. The lack of observable incommensurability 
for underdoped YBC06.33 (this work) and YBC06.35— 
is one of the differences that we find from the YBC06+X 
(highly doped) and LSCO families. This may indicate 
that the presence of stripes is not crucial for superconduc¬ 
tivity, although a breaking up of collinear spin order cer¬ 
tainly is. It would have been hard to detect incommen¬ 
surate correlations in a sample that is twinned and with 
disordered chains because the broadening would have ob¬ 
scured any anisotropy. Thus no attempt to do indepen¬ 
dent H and K scans was made to detect the cross-shaped 
scattering discussed by Gaw et af®. 

When comparing the properties of LSCO with YBCO, 
one must note the presence of intrinsic disorder caused by 
doping in LSCO®. The dopants occupy planes in close 
proximity to the active Cu0 2 planes but not in YBCO. 
Perhaps the ideal system for investigating the intrinsic 
properties of low-doped superconductors with minimal 
disorder effects would be cuprates with a larger num¬ 
ber of Cu0 2 planes 81 since doping occurs well away from 
most of these planes. In fact recently it was reported 
that in a Hg-based superconductor with five Cu0 2 lay¬ 
ers where the effects of disorder are minimal, long-range 
antiferromagnetism and superconductivity coexist!®. 

Our results are different from previous studies 12 — 
on detwinned low-doped YBC06.3 (T c =0), YBC06.35 
(T c =10 K) and YBC06.45 (T c =35 K) crystals where 
incommensurate elastic peaks with strong ab-plane 
anisotropy have been reported. These samples®*® were 
annealed for a short time (^1 day) and their SC transi¬ 
tion temperatures are lower than the UBC crystals stud¬ 
ied here which were annealed 18 20 for much longer (^3 
months). Different planar doping for similar oxygen con¬ 
tents, therefore, may be the cause of observed differences. 
For a similar linear dependence of the static incommen¬ 
surability vs. doping in LSCO along [1 0 0)-— -1®, the in¬ 
commensurability for YBCO with a doping of 0.055 (our 
sample) is expected to be 0.055 r.l.u. in H along [1 1 
0]. From elastic H-scans along [1 1 0] (Fig.|4|, we put an 
upper limit of 0.024 r.l.u. for possible static incommensu¬ 
rability in YBC06.33, certainly lower than the expected 
value. For the inelastic peaks, we note that the FWHM 
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FIG. 18. Magnetic excitations in undoped insulator (left), lightly doped SC YBC06.33 (middle) and optimally doped YBCO 
(right) are schematically depicted.). Spin-wave spectrum observed for the insulator can be described 4 by a Heisenberg Hamil¬ 
tonian with strong anisotropic superexchange interactions. 


of 0.058 r.l.u. observed in scans in H along [1 1 0] at 3 
meV, if it arose from incommensurate peaks along pure 
[1 0 0] direction, would correspond to 0.058 r.l.u. sep¬ 
aration of the putative incommensurate peaks. This is 
larger than the reported 13 value of 0.025 r.l.u. for the 
dynamic incommensurability at 3 meV in YBCO with 
similar oxygen content. 

Even if our inelastic data may agree with the pres¬ 
ence of dynamic incommensurability, from our elastic 
data, we can rule out the static stripe scenario for the 
very low doping of YBC06.33 (T c =8.4 K). This is in 
agreement with previous conclusion^ 22 2 ^ for the slightly 
higher doped YBC06.35 (T c =18.5 K). It appears that 
for very low concentrations of hole doping in YBCO, an 
arrangement of antiferromagnetically correlated quasi¬ 
static spins in three dimensions over a finite range is more 
favored than alternating quasi-one-dimensional spin and 
charge regions. This may be due to the suppression of 
superconductivity compared to magnetism or the pres¬ 
ence of only short-range oxygen chain order in very low 
doped YBCO. 

We find that the YBC06.33 spin dynamics are re- 
laxational with a relaxation rate, T, that saturates at 
the lowest temperatures. At high energies the dynamic 
planar correlations lengths shorten with an energy to q- 
width ratio similar to the spin-wave velocity in the insu¬ 


lating parent compound. Such a response in the energy 
domain (Eq. [t|) together with the Lorentzian momentum 
dependence (Eq. [5| describes the local relaxation of an¬ 
tiferromagnetically coupled spins whose memory is lost 
after a characteristic correlation time (1/T) and over cor¬ 
responding correlation lengths in the ab-plane and along 
the c direction. The finite relaxation rate as T—>>0 reflects 
how the doped system avoids Neel condensation by hav¬ 
ing spin waves scatter off locally ferromagnetic regions 
surrounding doped holes. Therefore, this low energy sup¬ 
pression may be regarded as necessary for removal of spin 
fluctuations that might compete with superconducting 
order. There are thereby fewer pair-breaking spin fluctu¬ 
ations that match superconducting gap energies of order 
3.5 k B T c =2.5 meV. 

The analysis of the temperature dependence of the 
magnetic excitations and the central mode indicates that 
^30 K can be regarded as the temperature scale for the 
appearance of the 3D spin-glass behaviour in YBC06.33. 
Below this temperature not only does the strength of 
the damped fluctuations deviate from a Curie-Weiss law 
(Fig.[8|b)) but also the central mode becomes more static 
and three-dimensional (Fig. 0d,e ,f)). This temperature 
is also close to the temperature below which quasistatic 
correlations gradually grow on cooling when measured 
with the very high energy resolution of backscattering 
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(see Fig.^b)). 

Even though we observe the development of magnetic 
correlations at the AF zone centre on cooling, their 
damped spin response stands in stark contrast to the 
low-energy spin waves in the insulator and to the hour¬ 
glass dispersion of more highly doped superconductors 
(see Fig. 18 for a schematic comparison). While doping 
appears to have only a modest effect on the high energy 
excitations. The insulator undergoes a Neel transition to 
3D long-range AF order where both the elastic scattering 
and the spin-waves can be described^ by an anisotropic 
exchange Heisenberg Hamiltonian. For our system, we 
observe instead a gradual increase of static short-range 
magnetic correlations with no low-energy spin waves. For 
higher doped YBC06+X with x>0.45, the doped holes 
frustrate any AF order causing the entire scattering to 
be dynamic (no elastic magnetic scattering). The spec¬ 
trum consists of a prominent resolution limited peak at 
the resonance energy, E res , locate d at the saddle point 
of the hourglass-shaped dispersioiP®^®^. The ob¬ 
served hourglass magnetic spectrum and presence of in¬ 
commensurate peaks around AF zone centre in cuprates 
are explained in terms of collective modes of an underly¬ 
ing stripe sta t<PEZI 

Since the only feature observed is the broad scatter¬ 
ing at low energies (~3 meV), one might conclude that 
this overdamped Lorentzian is the low-energy remnant 
of a resonance in YBC06.33. In fact from the empiri¬ 
cal relatiorP^ of E res =5.3k£T c and a T c of 8.4 K, one 
might expect the resonance excitation to occur at 3.9 
meV for YBC06.33. Moreover, providing the integrated 
intensity of the resonance is related to the superfluid 
density^, one might argue that because YBC06.33 is 
such a weak superconductor (T c ^8.4 K), the intensity 
of the resonance may be too small to observe. The mod¬ 
ified Lorentzian we observe, however, does not display 
the properties typical of a resonance!^ 20 KMESI neither an 
increase in its intensity below T c , nor an hour-glass dis¬ 
persion nor a decrease with applied magnetic field (see 
Section HID). Instead the overdamped exciations decline 
on cooling as the spins drop into the static ground state 
revealed by the growing central peak. We, therefore, re¬ 
gard the overdamped scattering peaked at ^3 meV as 
the spectral form of the excitations of an organized spin 
ground state with short-range correlations and incipient 
3D response. Similar spin-glass behaviour is observed in 
lightly doped La2Cuo.94Lio.o604 5 ^, Lai.gsBao.osCuCH 54 
and La2Cuo.95Zno. osOjsa 


A key point is that neither the growth in correlation 
length, the intensity of the central mode nor the spin 
dynamics show an anomaly at the superconducting tran¬ 
sition temperature T c =8.4 K, thus suggesting the spin 
correlations ignore superconductivity. This is fundamen¬ 
tally different from the higher doping (x>0.4-0.5) region 
of the phase diagram. Since the transition to the super¬ 
conducting state has no effect on the spin response, it 


appears that superconducting pairing of charges evolves 
independently of the glassy spin state as 3D magnetic 
correlations begin to dominate. This could be due to the 
fact that the pairs are already formed at temperatures 
much higher than T c or that the electrons responsible 
for superconductivity are different from those responsi¬ 
ble for the magnetic scattering. The site-based Neel spin 
pairing evolves as the momentum-based spin pairing of 
charges is developing. 

The stark difference we see between the properties 
of the low-doped yet superconducting YBC06.33 and 
the higher doped superconductors may be related to the 
presence of a critical hole concentration for a metal-to- 
insulator crossover (MIC) in the superconducting dome 
as suggested earlier 19 from neutron scattering studies 19 
on YBC06.45. The magnetic excitations were found to 
be gapless with a resonance much broader in energy than 
at the higher doping. A fundamental change in electronic 
properties within the SC dome can be similarly inferred 
from NMRf°°l and transport measurements in a field 1 () 1 1 
The c-axis transport measurements^ have also indicated 
that the coherence temperature for c-axis conductivity 
extrapolates to zero for doping close to p ~ 0.08. Al¬ 
though charge order has now been observed^] above a 
10 T critical field in the vicinity of x=0.5 in YBC06+X 
and fo r p= 0.11 where the Ortho-II chain structure is 
presentP^, there is no X-ray report for charge order at 
the low oxygen contents such as in the YBC06.33 we 
have studied. 

Although observations are qualitatively similar to 
higher doped YBC06.35 (T c =18 these are the 

first measurements on YBC06.33 with such a low T c 
revealing that even for such low doping, long-range anti¬ 
ferromagnetic order and superconductivity do not coex¬ 
ist. Rather the state coexisting with superconductivity is 
best regarded as a textured spin phase with short-range 
AF correlations. The measurements presented here on 
YBC06.33 also enable us to predict that a transition to 
a true 3D AF state will most likely occur only outside 
the dome at a doping below p c . 

Our results demonstrate that spin properties of very 
low-doped superconductors are fundamentally different 
from both the higher doped superconductors and the un¬ 
doped insulators. Any microscopic theory for HTSC su¬ 
perconductivity must account for the strong doping de¬ 
pendence of the observed properties of spin dynamics. 
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APPENDIX: ABSOLUTE INTENSITY SCALE 


The measured magnetic scattering depends on the sample size and can be put on absolute scale using the observed 
integrated intensity of an acoustic phonon. This is because the cross section of acoustic phonon scattering in the 
long wavelength limit is known and can be easily calculated. In this experiment we measured the transverse acoustic 


phonon close to (006) at (-0.15 -0.15 6) at 85 K as shown in Fig. 19 



FIG. 19. The transverse acoustic phonon measured at (-0.15 -0.15 6) close to nuclear Bragg peak (006) at 85 K. The fit is to 
a Gaussian function. The data shown is corrected for the presence of higher harmonics in monitor count rate. The integrated 
intensity of the peak above background is 3.68xl0 -7 (eV counts/mon). 


The energy-integrated intensity, J p h(Q)= f duo /(Q, cj), for a coherent one-phonon creation process can be written 104 


as: 


WQ) =AN ICv(Q)| 2 e “ 2W/(Q) cos 2 /? [1 + n(w p )] (8) 

where A is the scale factor, N is the number of nuclear unit cells, huo p is the phonon energy, [1 + n(uo p )] = [1 — 
e -hu>/k B T] 1 j g g ose f ac tor, e~ 2W ^ is the Debye-Waller factor (=1 in our calculations), M is the mass of the 
unit cell, |Fjv(Q)| is the static nuclear structure amplitude of the Bragg reflection nearest to where the phonon is 
measured, and (3 ^0 is the angle between Q and the phonon eigenvector. From the data presented in Fig. [19] and 
Eq-i we find A x A=2.04xl0 18 (cm 2 /eV-counts per mon). The constants used to obtain this ratio are given in 
Table QH 

The observed magnetic scattering intensity, / m ag(Q?^) is directly proportional to the magnetic scattering cross- 
section (see Eq. [l]) and can be written^! as: 

WQ,w) = A N mag g 2 (f) 2 |/(Q)| 2 B 2 (Q c ) e~ 2W S( Q,w) (9) 

where A is the scale factor, 7V mag is the number of magnetic unit cells (=N n so that the absolute intensities are deter¬ 
mined per formula unit), g ^ 2 is the Lande factor, 7 n =1.913 the gyromagnetic ratio of the neutron, ro=2.817xlO -13 
cm the classical electron radius, e~ 2W ^ is the Debye-Waller factor (=1 in our calculations), |/(Q)| is the Cu 2+ 
anisotropic magnetic form factor^, and B(Q C ) is the bilayer structure factor in YBCO giveiF 3 by 

B(Qc) = 2 sin(^^ u2 ) (10) 

where Q c = 2ixL/c and z' Cu2 = 1 — 2zq U 2 is the intra-bilayer spacing. The scattering function, 5 (Q,cj), is the Fourier 
transform of the spin-spin correlation function and is related to the spin susceptibility through the fluctuation- 
dissipation theorem 
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TABLE II. The constants use d in the phonon absolute scale calculations. The transverse acoustic phonon data used in the 


calculation is shown in Fig. 16 As discussed in the text cos 2 /3=l. We have also approximated e 2W ^ to 1. Lattice parameters 


used are a=3.844 A, b=3.870 A, and c=11.791 A 


Q 

Q 2 

|Fiv(Q)| 2 

M 

ftw p 

1 + n{uj p ) 

4°t s> (Q) 

(r.l.u.) 

(cm 2 ) 

(cm 2 ) 

(kg) 

(meV) 

(unity) 

(eV counts/mon) 

(0 0 6) 

I.02xI0 17 

I.73xl0“ 23 

1.09xl0“ 24 

5.76 

1.84 

3.68xl0 -7 


S( Q,w) 


[1 + «M1E 


a, /3 





( 11 ) 


with the summation over Cartesian directions. Since there is no evidence for the presence of a long-range magnetic 
order and a preferred orientation of the moments, the summation in Eq. [TT] must be invariant with respect to the 
rotation of the indices. Therefore, the sum will be equal to (2x"/7t (gfis) 2 ) where x" is the isotropic susceptibility per 
formula unit. Using this result, Eq. [9] can be rewritten as 


WQ>") = A N mag g 2 (Y) l/(Q)| 2 


B 2 (Qc 


,-2W(Q) 


[1 + n(o 


2x"(Q^) 

ir{gLiB ) 2 


( 12 ) 


The constants used in the calculation of the absolute intensity at (0.5 0.5 2) are lis ted i n Table |HI| The observed inten¬ 
sities of constant-Q energy scans, can now be put on an absolute scale using Table III and A x A’=2.04xl0 18 (cm 2 /eV- 
counts per mon) obtained above 


[1 + n(w)] x"(Q, w) = (3.38 x 10 6 ) x / mag (Q, w) (13) 

where x"(Q,w) is in (/j|,/eV) per unit-formula and J mag is in counts/mon. The absolute scales in Fig. [sja) are 
obtained from this method. 


TABLE III. The constants used in the calculation of the absolute intensity of the magnetic scattering measured at (0.5 0.5 2). 
The number of magnetic unit cells is set equal to the nuclear unit cells N mag =N so that the absolute intensities are determined 
per unit formula. 


Q 

z Cu2 

(W 

l/(Q)l 2 

B\Q C ) 

(r.l.u.) 

(unity) 

(cm 2 ) 

(unity) 

(unity) 

(0.5 0.5 2) 

0.36 

7.27xl0 2e 

0.81 

3.86 
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